Generic mixed columnar-plaquette phases in Rokhsar-Kivelson models 
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We revisit the phase diagram of Rokhsar-Kivelson models, which are used in fields such as su- 
perconductivity, frustrated magnetism, cold bosons, and the physics of Josephson junction arrays. 
From an extended height effective theory, we show that one of two simple generic phase diagrams 
contains a mixed phase that interpolates continuously between columnar and plaquette states. For 
the square lattice quantum dimer model we present evidence from exact diagonalization and Green's 
function Monte Carlo techniques that this scenario is realised, by combining an analysis of the ex- 
citation gaps of different symmetry sectors with information on plaquette structure factors. This 
presents a natural framework for resolving the disagreement between previous studies. 



The Rokhsar-Kivelson Quantum Dimer model (RK- 
QDM) m on the square lattice, originally proposed in the 
context of high-temperature superconductivity, and its 
descendants have taken on a central role in the study of 
quantum systems incorporating a hard local constraint. 
They have thus been prominent in the context of hard- 
core bosons hopping [3| on frustrated lattices, Josephson 
junction arrays 0], frustrated Ising models in a trans- 
verse field or with small XY exchange [3| , gauge theories 
in unusual sectors [1], spin orbital models U and cold 
atoms [?[• Studies of this model and its extensions have 
unearthed a wealth of phenomena, including instances 
of deconfincd quantum criticality and a new route to 
deconfinement [§] . 

As many RK models belong to the rare class of mod- 
els of correlated quantum matter without a sign prob- 
lem, they are in principle amenable to efficient numerical 
study; they also tend to have well-studied effective field 
theories formulated in terms of height/gauge degrees of 
freedom This makes it all the more surprising that 
there is still substantial disagreement about the phase 
structure of the original RK model, the QDM on the 
square lattice. Indeed, a pioneering study by Leung et 
al.[l3| suggested a transition out of a columnar phase 
(Fig.[T]) to occur as v/t increases from — oo at v/t ^ —0.2 
(see Eq. [1]); however, a careful detailed recent investiga- 
tion by SyljuasenQ argued that in fact the columnar 
phase persists until v/t '-^ -1-0.6 

A direct columnar-plaquette transition en route to the 
RK point at V = t IS in fact not unusual - it appears to 
exist in the closely related quantum six-vertex model on 
the square [l^l and the hexagonal QDM 11 1. Here, we 
show that there exists a second generic phase diagram 
where this first-order transition is replaced by a contin- 
uous interpolation via a mixed phase (Fig. [T|). Based 
on Exact Diagonalisations (ED) and Green's Function 
Quantum Monte Carlo (GFMC) numerics we identify 
the square lattice QDM as the first candidate for real- 
ising this scenario: around w/t = 0, we find evidence of 
a mixed columnar-plaquette phase, which breaks trans- 
lational symmetry in two perpendicular directions like 



the plaquette, but also 7r/2 rotational symmetry like the 
columnar phase. The strongest evidence for such a phase 
is provided by a symmetry-based finite-size scaling anal- 
ysis of the low energy spectrum. 

The remainder of this paper provides the analysis back- 
ing up these assertions: first we present the pertinent nu- 
merical results on the QDM, followed by a more general 
analysis of RK models in d = 2 whose effective theory 
can be described by a height model. 

The square lattice QDM Hamiltonian reads: [ij 



H = vY,Nf{c)\c){c\~tY,\c){c'\ 



(1) 



(c,c') 



where the sum over c runs on all configurations in the 
Hilbert space, Nf(c) is the number of flippable plaque- 
ttes contained in |c), i.e. number of plaquettes with two 
parallel dimers, and the sum over (c, c') runs on all con- 
figurations |c) and |c') that differ by a single dimer flip. 
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FIG. 1: Schematic representation of the different or- 
dered states considered, the columnar, the mixed columnar- 
plaquette and the plaquette states. Dashed lines show the ex- 
tra bonds of the checkerboard lattice of the equivalent hard- 
core boson model (see Ref. 01). On the left, dimers (thick 
bonds) map to bosons (red circles). 

Symmetry and first excitation spectrum: To compute 
excitation gaps we combine ED with GFMC calculations 
of dynamical correlations in imaginary time. In the lat- 
ter case, the choice of an operator is crucial in order 
to "construct" an excitation with well-defined symme- 
try properties (i.e. quantum numbers). Starting with the 
elementary symmetry- breaking patterns of Fig. [H degen- 
erate ground states with different quantum numbers can 
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be constructed as listed in Table |T1 On general grounds, 
symmetry breaking involving any of these states is sig- 
naled on finite clusters of increasing size by the collapse 
of the quasi-degenerate GS and the opening of a robust 
gap above them to higher energy excitations. Clearly 
columnar and plaquette phases differ in two states with 
different quantum numbers (defined by a momentum and 
a point group quantum number), (F, Bi) (point q = (0, 0) 
in the BZ depicted in the inset of Fig. [3l with symme- 
try Bi) for the first one and {K,Ai) {q = (7r,7r)) for 
the second. As shown in [l^], a level crossing between 
two such energy levels was the signature of a first order 
phase transition. In contrast, as we show further down, 
no such level crossing is observed here, suggesting a more 
complex scenario. 
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TABLE I: Quantum numbers of the degenerate GS charac- 
terizing the ordered phases considered in this paper. We used 
the standard notation of irreducible representation of the Cav 
and C2v point group symmetries, acting around the site O as 
depicted in Fig.[T] Definition of the momentum points is given 
in the Brillouin zone of the inset of Fig. [S] Note that a state 
at momentum (0, tt), not shown, is degenerate with (tt, 0). (*) 
denotes the first excited level in the {M,Ai) sector. 

Since ED are limited to 8 x 8 clusters we also used 
GFMC on clusters of size up to 22 x 22 to calculate gaps 
in different symmetry sectors by considering dynamical 
correlations as defined in lisll: 



Diq, 



{^G\Pc.i-q)e-"^Pa{q)\^o) 
(*G|e-«1*o) 



(2) 



with I^'g) a guiding function. To improve the quality of 
5*0) we explicitly work in the boson representation (see 
Fig. [T] and Ref. [3] ) which enables to minimize a Jastrow 
wave-function with the help of a variational Monte-Carlo 
method. Pa{q) is a (diagonal) operator with the same 
symmetry as the excited state that we intend to target 
and is defined as the Fourier transform of plaquette op- 
erators 



k ■' I 



= d{i)d{j) + d{k)d{l) 
= d(i)d{j) - d{k)d{l) 



(3) 



where d{i) is the dimer operator acting on link i. Note 
that the point group symmetry of Pa depends on the 
momentum q. For P+, it always corresponds to the most 
symmetric Ai irreducible representation (IR) . For P_ , it 



corresponds to the Bi IR for the high-symmetry points, 
e.g. r and and to the Ai IR at point M . Note also 
that the expectation value of P^{K) (P_(r)) is finite in 
a generic plaquette (columnar) state. 
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FIG. 2: First gap excitations, a) Comparison between ED and 
GFMC for the 8x8 site cluster for the interesting q-points. 
Size-scaling of the gaps for the {K,Ai) (b) and (F, Bi) (c) 
symmetries. 

Our results for the excitation gaps are given in Fig. [2l 
with q-points defined through the Brillouin zone dis- 
played in Fig. [3] Panel (a) shows the full exact spec- 
trum on a 8x8-size cluster and a comparison with the 
GFMC estimates for the lowest levels. We find an ex- 
cellent agreement which validates the GFMC method on 
larger clusters. Interestingly, our data distinguish three 
regimes: (i) for v/t < -0.5, the (F,Si) and (M,Ai) 
levels have very small excitation energy and are well- 
separated from the rest of the spectrum by a sizable gap; 
(ii) for roughly v/t > —0.5, a second group of states 
of {K,Ai), {M,Ai) and {K,Bi) quantum numbers joins 
the previous low energy states while remaining well sep- 
arated from other higher energy states; (iii) the critical 
regime near the RK point u < i. A comparison with Ta- 
ble [T] immediately suggests a mixed phase and a colum- 
nar phase, respectively, for v/t > —0.5 and v/t < —0.5. 
These data seem already to be inconsistent with a pure 
plaquette state for which the ABi(r) (AAi(i^)) gap is 
finite (vanishing), which would imply a crossing of the 
(F, Bi) and {K^Ai) levels, in contradiction to what is 
seen numerically. To back up the above, we have per- 
formed finite size scalings of A^^ {K) and A^^ (F) shown 
in Fig.^b) and (c), respectively. A^^ {K) clearly extrap- 
olates to zero in the thermodynamic limit for v/t > 0.1 
(compatible with both plaquette and mixed states) and 
opens up near v/t ~ 0.0 ±0.1. Surprisingly, Asi(r) van- 
ishes in the thermodynamic limit in the whole range of 
parameters studied (within an accuracy of lO"'^) with a 
clear exponential behavior for v/t < —0.4. Note that the 
fact that ABi(r) always remains smaller than Aai{K) 
for all sizes itself rules out a pure plaquette state. 
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Structure factors: Next we check that the above sce- 
nario obtained from the spectrum analysis (beheved to 
be the most accurate criterium) is compatible with the 
behavior of the related structure factors. 



(4) 



(*ol*o) 

Results are displayed in Fig. [3] where the left panel shows 
the behaviors of the P+ and P_ correlations along a given 
path of the Brillouin zone (on a 8x8-size cluster) and the 
right panel the sizc-scalings of their respective order pa- 
rameters M{q) = yJl(ci)lL, L being the linear size of the 
system. These data correspond tov/t = 0.2, i.e. far from 
the proposed plaquette-columnar transition in Q , deep in 
the previously-supposed columnar phase. However, our 
results reveal a Bragg peak at point K which is found to 
survive upon extrapolation to the thermodynamic limit. 
On the right panel of Fig. [3l size-scalings suggest that 
columnar order only develops for v /t ^ 0.6 while a finite 
plaquette order is present in the range 0.0 < v/t < 0.8. 
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FIG. 3: a) Behaviour of the structure factor for the 8x8-site 
cluster at v/t = 0.2, along the path of the Brillouin zone 
depicted in the inset. A clear peak appears at point K. Size- 
scaling of the order parameters in the two symmetry sectors 
are given, for {K,A-i_) (b) and {T,Bi) (c). 

Discussion and interpretation of the numerics: Our 
extrapolated results as a function oi v/t for the rele- 
vant gap and structure factors are summarized in Fig. |31 
The results for the gaps provide clear evidence for 
lattice-rotation symmetry breaking for all parameters 
(ABi(r) = 0) up to v/t = 0.6 (for v/t > 0.6 a tiny (r,Si) 
gap cannot be excluded). We find two regimes: (i) For 
v/t < 0.0, Aai{K) is finite, indicating a pure columnar 
phase; (ii) for 0.0 < v/t < 1, Aai{K) vanishes, implying 
an additional translation-symmetry breaking compatible 
with the mixed phase. The behavior of the structure 
factors is consistent with such a picture in this range of 
parameters. However, (i) error bars for —0.1 < v/t < 0.1 
are still too large for simultaneous accurate extrapola- 
tions of both gaps and order parameters and (ii) upon 




FIG. 4: Thermoynamic limit of the order parameters (black 
and red curves) and of the first excitation gaps (blue and green 
dashed curves) in function of v/t. Mai (K) is multiplied by a 
factor 5 for clarity. A zoom around v/t = 0.6 is shown in the 
inset. Mai{K) is not multiplied here. 



approaching the continuous transition at the RK point, 
both structure factors become so small (inset in Fig.d]), 
that we are unable to distinguish them from zero within 
the accuracy of our simulations for v/t ^ 0.8. 

Syljuasen argued that rotation symmetry is restored 
when v/t > 0.6 signaling the transition to a pure pla- 
quette phase, at least up to v/t ~ 0.9 This claim 
is consistent with our results in that a small (T,Bi) gap 
cannot be ruled out for v/t > 0.6 and that it is plausible 
that our diagonal plaquette operator, although of the 
correct symmetry to "pick up" plaquette order, only gives 
a very small weight to the dynamical correlations when 
approaching the RK point. Also, note that the transi- 
tion at v/t = 0.0 ± 0.1 that we evidence here from the 
columnar state to the novel mixed columnar-plaquette 
state seems to be smooth and might be compatible with 
a second order phase transition. We believe our scenario 
of an intermediate mixed state reconciles the a priori 
conflicting results, (i) the ED data by Leung et al. sug- 
gesting the existence of an intermediate phase for roughly 
v/t > —0.2 [la] and (ii) Syljuasen's claim of a transition 
at v/t = 0.6 [9|]. However, Syljuasen's analysis suggested 
that a mixed phase, if present at all, exists in a narrower 
range of parameters. At the same time, he pointed out 
that there was no sign of a discontinuous transition. 

Effective height theory: We next show that the generic 
alternative to a first-order plaquette-columnar transition 
is a continuous interpolation between them via a mixed 
phase, in agreement with the scenario obtained above. 
This we do by reformulating the dimer model in terms 
of a height variable, living on the plaquettes of the 
square lattice. Mutatis mutandis, such a mapping applies 
to d = 2 RK models incorporating constraints which can 



be cast as local U(l) conservation laws |10l. Illl Il4 



The mapping from dimcrs to heights proceeds as fol- 
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lows. Label the two sublattices of the square lattice as 
red and green. When going from the centre of one pla- 
quette to another in an (anti-)clockwise direction around 
a (red) green site, (add) subtract 1 from the value of the 
height if no dimer is crossed; and (subtract) add 3 if a 
dimcr is crossed. An effective field theory starting from 
the microscopic heights can then obtained for a coarse- 
grained height field h. This procedure is described in 
detail in a nice paper by Zeng and Henley [3], and the 
corresponding quantum theory is given in Ref. |5|. 

For our purposes, the following facts are important. 
Firstly, a 'flat surface', i.e. long-range order in h, corre- 
sponds to a dimer crystal. In particular, a (half-) integer 
valued {h) corresponds to a (columnar) plaquette phase, 
whereas intermediate values of (/i) correspond to mixed 
phases like the one discussed above. Secondly, h h+1 
amounts to a 7r/2 rotation of the dimer configuration, so 
that h ^ h + 4 leaves the dimer configuration unchanged. 
Finally, the leading terms in an effective action incorpo- 
rating all terms not ruled out by symmetry or microscopic 
considerations of the model in d = 2 + 1 dimension are 



S = 



d^xdr 



[{drhf + P2{\7hf 
+ A cos(27r/i) 4- A* cos(47r/i)] , (5) 



where r denotes Euclidean time. This action differs from 
that in Ref. H in two ways: (i) a term (V^/i)^ is missing 
because it plays an important role only at the RK point 
where p2 vanishes, whereas we are interested in what hap- 
pens elsewhere; (ii) the term p, cos(47r/i) has been added 
as its presence is necessary for the analysis of the transi- 
tion out of the columnar phase. To see this, let us min- 
imise S for uniform configurations. For /i = 0, the sign 
of A determines if one has a columnar (A > 0) or pla- 
quette state. At the transition point between the two, 
A = 0, implying a continuous degeneracy ((/i) can take 
any value £ [0,4[), which has so far not been observed. 

For fi 0, the fate of this transition depends on the 
sign of fi. For < 0, the minima of the A term coincide 
with those of the /x term, independently of the sign of A. 
This destroys the continuous degeneracy at A = 0, and 
leads to a simple first order transition between columnar 
and plaquette states. By contrast, for /x > 0, the and 
A terms compete: the minima of the latter coincide with 
the maxima of the former. Up to symmetries, one obtains 
(h) changing continuously between at A = — 4/Lt and 1/2 
at A = 4/x as 



, if A < -4^, 
27r(/i) = arccos(-A/4/i) , if |A| < 4^, 
TT , if A > 4^ . 



(6) 



1 before v/t does, the mixed phase terminates before 
the RK point; otherwise, it will be terminated by the 
RK transition. We note that other 'non-generic' scenar- 
ios (e.g. for particularly small values of y^t, or for non- 
monotonic dependence of A//x on v/t) are also possible. 

Concluding remarks: The analysis of low energy spec- 
tra offers a powerful method to investigate spontaneous 
symmetry breaking, in particular eliminating the bias 
that the choice of an order parameter entails. While our 
system sizes - up to 30 x 30 sites with GFMC - are much 
bigger than those usually accessible to ED, they are still 
subject to the uncertainties of extrapolating delicate or- 
dering phenomena to the thermodynamic limit. 

However, our analytics support the numerics by show- 
ing that the phase diagram incorporating an intermedi- 
ate mixed phase is one of two generic possibilities for such 
RK models, the other being a direct first order columnar- 
plaquette transition. The square lattice QDM is the first 
strong candidate for realising this new scenario. 
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In the region |A| < Afx we therefore find a low-symmetry 
mixed state which continuously crosses over from primar- 
ily columnar to plaquette character. If — A/4/i reaches 
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